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Beschreibung 

Verfahren zur Oberf lachenkonturierung eines dreidimensionalen 
Abbildes 

5 

Die Erfindung betrifft ein Verfahren zur Oberf lachenkonturie- 
rung eines dreidimensionalen Abbildes eines Objekts. 

Fur verschiedene Anwendungen, beispielsweise in der Pho- 
10 togrammetrie, in der industriellen Bildverarbeitung, bei der 
Uberwachung, im Verkehrswesen, in der Medientechnik und der 
^nimation, sind Bildverarbeitung und Obj ekter kennung zweidi- 
mensionaler und teilweise auch dreidimensionaler Abbilder von 
Objekten relativ weit vorangeschritten. Die die Abbilder um- 
15 fassenden Bilddatensat ze sind in der Regel relativ schwach 
gestort und sind z.B. mit optischen Sensoren aufgenommen 
(vgl. z.B. Inspect. MessTecAutomation, Sonderpublikation 
10/01, GIT-Verlag Darmstadt, 2001 oder Vision Systems Design, 
Vol. 6 #9, Pennwell, Sept. 2001). 

20 

Sind die Bilddatensat ze relativ schwach gestort, konnen die 
Oberf lachenkonturen der Abbilder der Objekte mittels bekann- 
ter Verfahren relativ einfach ermittelt werden. Mit linearer 
oder nichtlinearer Filterung kann ein eventuelles Rauschen 
: |0£5 der Bilddatensat ze weiter verringert werden. Urn Texturen der 
Abbilder zu erkennen, kann anschlieftend ein so genanntes Pat- 
tern Matching durchgefuhrt werden. Grenzen der Abbilder kon- 
nen beispielsweise durch 2D- oder 3D-Gradientenbildung ermit- 
telt werden. Auf Gradientenbildung basierende Techniken sind 

30 z.B. "active contours", beschrieben in Chalana V., Linker D., 
"A Multiple Active Contour Model for Cardiac Boundary Detec- 
tion on Echocardiographic Sequences", IEEE Trans. Med. Ima- 
ging 15, #3, Juni 1996, oder so genannte "snakes", beschrie- 
ben in Mclnerney T, Terzopoulos D: "T-snakes: Topology Adap- 

35 tive Snakes", Medical Image Analysis 4, Seiten 73-91, 2000. 
Abgebildete Objekte gleicher Helligkeit, Farbe oder Textur 
konnen andererseits z.B. durch das so genannte region growing 
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"ausgef ullt " und damit fur eine bestimmte Auf gabenstellung 
vollstandig beschrieben werden (vgl. Jendrysiak U., "Segmen- 
tierung von Schnittbildern" , Spektrum der Wissenschaft, Dos- 
sier 1/1999: Perspektiven in der Medizintechnik 24-29, 1999) 

5 

In der Medizintechnik kann es jedoch sein, dass abgebildetes 
gesundes und krankes Gewebe einen relativ geringen Kontrast- 
unterschied aufweist und Bilddatensat ze relativ stark ver- 
rauscht sind. Das gilt insbesondere fur Rontgenbilder (Compu 
10 tertomographie, Angiographies Fluoroskopie ) , mit einem Mag- 
netresonanzgerat oder in der Nuklearmedizin hergestellte Bil 
N^ler und insbesondere auf Ultraschall basierte Bilder (vgl. 
Marais P., Brady J. , "Detecting the Brain Surface in Sparse 
MRI Using Boundary Models", Medical Image Analysis 4, 283- 
302, 2000, Chalana V., Linker D. , "A Multiple Active Contour 
Model for Cardiac Boundary Detection on Echocardiographic Se 
quences", IEEE Trans. Med. Imaging 15, #3, Juni 1996, oder 
Sakas G., "Dreidimensionale Bildrekonstruktion aus Ultra- 
schall-Daten" , Spektrum der Wissenschaf t , Dossier 1/1999, 
Perspektiven in der Medizintechnik 18-24, 1999) . 



Insbesondere bei einer Diagnosef indung fur einen Patienten 
muss in der gesamten Tiefe des Korpers oder der Organe des 
Patienten nach feinen Ob j ektunterschieden dif f erenziert wer- 

^5'^ den. Da die mit medizinischen Geraten auf genommenen Bildda- 
tensatze oft relativ stark verrauscht sind und aufgenommene 
Strukturen relativ komplex sein konnen, ist bereits ein Er- 
mitteln einer Kontur in einem zweidimensionalen Bild relativ 
rechenauf wandig . Das Ermitteln einer Kontur eines dreidimen- 

30 sionalen Abbildes eines Objektes ist sogar noch schwieriger. 
Bekannte Verfahren wurden deshalb verfeinert und oftmals ad- 
hoc an die vorliegende Problemstellung angepasst. Ein Bei- 
spiel ist in Chalana V., Linker D., "A Multiple Active Con- 
tour Model for Cardiac Boundary Detection on Echocardi- 

35 ographic Sequences", IEEE Trans. Med. Imaging 15, #3, Juni 
1996 beschrieben. Relativ rechenauf wandig ist auch ein in 
Jendrysiak U., "Segmentierung von Schnittbilder " , Spektrum 
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der Wissenschaft, Dossier 1/1999, Perspektiven in der Medi- 
zintechnik, Seiten 24-29, 1999 beschriebenes Verfahren. 

Die Aufgabe der Erfindung ist daher, ein Verfahren anzugeben, 
5 mit dem eine Oberf lachenkonturierung eines dreidimensionalen 
Abbildes eines Objekts mit geringerem Rechenauf wand durch- 
fuhrbar ist. 

Die Aufgabe der Erfindung wird gelost durch ein Verfahren zur 
10 Oberf lachenkonturierung eines dreidimensionalen Abbildes ei- 
nes Objekts, aufweisend folgende Verf ahrensschritte : 




- Legen eines ersten Sets von Ebenen durch das dreidimensio- 
nale Abbild, 

15 

- Ermitteln der Konturen, die in jedem der Ebenen des ersten 
Satzes von Ebenen abgebildet und der Oberf lachenkontur des 
dreidimensionalen Abbildes zugeordnet sind, und 

20 - Zusammenf iigen der in jeder Ebene des ersten Sets von Ebenen 
ermittelten Konturen zu einem der Oberf lachenkontur des 
dreidimensionalen Abbildes zugeordneten Oberf lachengitters . 

Erf indungsgemafi soli die Oberflache eines dreidimensionalen 
vjjpo Abbildes eines Objekts konturiert bzw. segmentiert werden. 

Das dreidimensionale Abbild ist in einem Volumendatensat z , 
also einem Datensatz von 3D-B16cken B(x,y,z) mit KxLxM Voxeln 
(volume elements) gespeichert. Das dreidimensionale Abbild 
kann aber auch in Form eines Datensatzes gespeichert sein, 
30 der K Schichten von Bildern mit je LxM Pixeln (picture ele- 
ments) umfasst. 

Dazu werden die Ebenen des ersten Sets von Ebenen durch das 
dreidimensionale Abbild gelegt. Anschliefiend wird fur jede 
35 einzelne dieser Ebenen die Kontur, die der Oberf lachenkontur 
des dreidimensionalen Abbildes zugeordnet ist, ermittelt. Das 
Ermitteln der den einzelnen Ebenen zugeordneten Konturen ent- 



200211387 



20 
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spricht dem Ermitteln von Konturen zweidimensionaler Abbil- 
der. Danach werden die einzelnen Konturen zu einem die Ober- 
f lachenkontur des dreidimensionalen Abbildes zugeordneten O- 
berf lachengitters zusammengef iigt . 

Nach einer bevorzugten Ausf uhrungsf orm der Erfindung schnei- 
den sich die Ebenen des ersten Sets von Ebenen in einer ers- 
ten Gerade, die gemaft einer besonders bevorzugten Variante 
der Erfindung wenigstens im Wesentlichen durch den Volumen- 
schwerpunkt des dreidimensionalen Abbildes verlauft. Dies ist 
besonders bei einem konvexen dreidimensionalen Abbild vor- 
/teilhaf t . 



Weist der das dreidimensionale Abbild umfassende Volumenda- 
tensatz beispielsweise einen relativ schwachen Kontrast auf, 
so kann insbesondere fur eine automatische Lokalisation des 
dreidimensionalen Objektes gemaft einer Variante der Erfindung 
fur das Ermitteln des Volumenschwerpunkts des dreidimensiona- 
len Abbildes der Volumendatensat z einer Filterung unterzogen 
werden. Die Filterung kann insbesondere eine Medianf ilterung 
sein und der Volumenschwerpunkt kann insbesondere automatisch 
ermittelt werden. Der Volumendatensat z wird dadurch unter 
Kantenerhaltung stark geglattet. Durch eine derartige starke 
Glattung verliert man zwar die Detaillierung der Begrenzun- 
gen, jedoch kann der Volumenschwerpunkt des dreidimensionalen 
Abbildes relativ klar bestimmt werden. 

Ist das abzubildende Objekt, wie es nach einer Variante der 
Erfindung vorgesehen ist, ein Teil eines Lebewesens und wird 
das dreidimensionale Abbild mittels eines medizintechnischen 
Gerates aufgenommen, so kann, wie es nach einer weiteren Aus- 
fiihrungsform der Erfindung vorgesehen ist, der Volumenschwer- 
punkt des dreidimensionalen Abbildes wahrend einer navigati- 
onsgefuhrten Behandlung des Lebewesens ermittelt werden. Der 
Volumenschwerpunkt kann beispielsweise wahrend der navigati- 
onsgefuhrten Behandlung durch die Position des Fokus eines 
navigierten medizinischen Instrumentes ermittelt werden. 
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Eine fur die Navigation geeignete medizinische Vorrichtung 
ist u.a. in der DE 199 51 502 Al beschrieben. Beispielsweise 
mit einem Rontgengerat wird ein Volumendatensat z von dem in- 
teressierenden Bereich des Korpers eines Patienten erstellt. 
5 Insbesondere wahrend eines minimal invasiven Eingriffs fiihrt 
z.B. ein Chirurg das medizinische Instrument in den Korper 
des Patienten ein. An dem medizinischen Instrument ist ein 
Positionssensor des Navigationssystems angeordnet, so dass 
das Navigationssystem die Position, d.h. die Lage und Orien- 

10 tierung des medizinischen Instruments, relativ zum Korper des 
Patienten bestimmt. Aufgrund der Positionsbestimmung, also 
^aufgrund einer Bestimmung der Ortskoordinaten des medizini- 
schen Instruments, kann anschlieftend ein Abbild des medizini- 
schen Instruments in das dem Volumendatensat z zugeordneten 

15 Bild, das z.B. mit einem Monitor dargestellt wird, eingeblen- 
det werden. 

Eine automatische Erkennung der ungefahren Position des medi- 
zinischen Instruments gelingt auch in der Nuklearmedizin . Ein 
20 mit einem PET auf genommener Volumendatensat z weist einen re- 
lativ guten Kontrast auf, hat aber eine relativ schlechte 
Auflosung. Eine weitere Alternative ist es, den Volumen- 
schwerpunkt manuell in einer anderen Darstellung zu bezeich- 
nen, z.B. im Rontgen unter Kontrastmittel in zwei Projektio- 
^5 nen oder in einer CT-MPR-Darstellung mit 1 Klick (MPR bedeu- 



tet "Multi-Planare Ref ormatierung" . Dies ist nichts anderes, 
als die Neuzusammenstellung des Volumendatensat zes in anderer 
Orientierung als z.B. ursprungliche horizontale Schichten 
(orthogonale MPRs : horizontal, sagittal, koronal; freie MPR: 

30 beliebig schrage Sekundarschichten) . Dies gelingt mit neuen 

Mehrzeilen— Spiral-CTs oder MR-Geraten isotrop z.B. aus 512 CT 
Schichten mit je 512 2 Pixeln) . Der Hintergrund dazu ist, dass 
verschiedene Ansichten des Abbildes des Objekts und seiner 
Umgebung gekoppelt in definierter Position und Orientierung 

35 dargestellt werden. Wird die dreidimensionale Ansicht ver- 
schoben und/oder gedreht, geht beispielsweise die entspre- 
chend dazu orientierte MPR-Darstellung nebenan mit. Die drei- 
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dimensionale Verschiebung und Drehung werden so lange fortge- 
flihrt, bis das interessierende Abbild von der MPR-Schicht 
zentral geschnitten wird. Dann wird in die Mitte dieser ebe- 
nen Darstellung des Abbildes geklickt; damit ist das Objekt- 
5 zentrum im Dreidimensionalen definiert. 

Mittels des erf inderischen Verfahrens wird die Dimensionali- 
ty der Volumenbearbeitung ( 3D-Bildbearbeitung) in wenigstens 
einer Stufe reduziert. Den Verlauf einer 3D-Ob j ektdichte D 

10 des dreidimensionalen Abbildes entlang einer Koordinatenachse 
x des Raumes kann man als zweidimensionale Funktion Dl (x) 

pf^N^graphisch darstellen. Das dreidimensionale Abbild im Raum 

kann daher dementsprechend mit einer Funktion D3(x,y,z) be- 
schrieben werden, deren Werte (Auspragungsvariable) in der 

15 vierten Dimension aufgetragen werden. Es wird nun hier ein 

Weg beschrieben, wie der von x, y und z aufgespannte Raum so 
projiziert werden kann, dass jeweils nur das Problem ahnlich 
der "eindimensionalen Irrfahrt" (vgl. Tittmann P., "Irrfahr- 
ten in Graphen", Memo Hochschule Mittweida, Sept. 1998.) bzw. 

20 eine Optimierung allein entlang einer Abszisse x' (statt (x, 
y, z) ) zu losen ist. 



Urn das zu erreichen, werden die Ebenen des ersten Sets von 
Ebenen durch das dreidimensionale Abbild gelegt. Die Ebenen 
•^p?-^ des ersten Sets von Ebenen sollten im allgemeinen Fall mog- 
lichst nicht redundant alle Voxel des dreidimensionalen Ab- 
bildes in verschiedenen Orientierungen, aber doch vollstandig 
uberdecken. Vorzugsweise schneiden sich die Ebenen des ersten 
Sets von Ebenen in der ersten Gerade, die wenigsten im We- 

30 sentlichen durch den Volumenschwerpunkt des dreidimensionalen 
Abbildes verlauf t. Die einzelnen Ebenen des ersten Sets von 
Ebenen konnen beispielsweise aquidistante Winkel oti bis a™, 
wenn das erste Set von Ebenen m Ebenen umfasst, bezogen auf 
eine feststehende Achsennormale, aufweisen. Durch den Uber- 

35 gang in eine Rotationsgeometrie ergeben sich Redundanzen in 
Richtung Drehpunkt und geringere Abtastung nach distal. Eine 
Abtastung braucht jedoch haupt sachlich im Bereich der erwar- 



200211387 



7 

teten Oberf lachenkontur beachtet zu werden. Durch die "Neuer- 
fassung" des dreidimensionalen Abbildes in gegeneinander ro- 
tierten Ebenen, die das dreidimensionale Abbild schneiden, 
erreicht man eine erste Stufe einer Dimensionserniedrigung, 
da die den Ebenen zugeordneten Datensatze zweidimensionale 
Datensatze sind und beispielsweise mit kartesischen Koordina- 
ten beschrieben werden. Wird z.B. der Volumendatensat z mit 
kartesischen Koordinaten (x,y,z) beschrieben, so kann jeder 
den Ebenen des ersten Sets von Ebenen zugeordnete Datensatz 
beispielsweise mit zweidimensionalen kartesischen Koordinaten 
(Xj*\yj* , J beschrieben werden, wobei j = l ... m ist, wenn das 
^erste Set von Ebenen m Ebenen umfasst. Anschliefiend werden 
die den einzelnen Ebenen des ersten Sets von Ebenen zugeord- 
neten Datensatze analysiert, urn Konturen, die in jedem der 
Ebenen des ersten Satzes von Ebenen abgebildet und der Ober- 
f lachenkontur des dreidimensionalen Abbildes zugeordnet sind, 
zu ermitteln. 

In jeder dieser Ebenen (Xj * \ Yj yy ) wird eine geschlossene 
Kontur erwartet. Die Kontur Kj w der j-ten Ebene des ersten 
Sets von Ebenen hat die Punkte K j± " = (Xj fi *\ yj,i*'), mit i 
= 1 ... n. Urn jedes Kji^ zu bestimmen, benotigt man also ei- 
ne zweidimensionale Suche in (Xj ' \ yj A M , die dem Fachmann 
im Allgemeinen bekannt ist. 

Urn den Rechenauf wand zu verringern, ist es jedoch gemaft einer 
besonders bevorzugten Ausf uhrungsf orm der Erfindung vorgese- 
hen, dass, wenn die Bildinhalte jeder Ebene des ersten Sets 
von Ebenen mit kartesischen Koordinaten beschrieben sind, fur 
das Ermitteln der Konturen jeder einzelnen Ebene folgende 
Verf ahrensschritte fur jede Ebene auszufiihren: 

- Durchfuhren einer Koordinatentransf ormation nach Polarkoor- 
dinaten, bezuglich des Volumenschwerpunktes des dreidimen- 
sionalen Abbildes und damit Abrollen der Kontur und 

- Ermitteln der Kontur in der transf ormierten Ebene. 
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Die Kontur in der transf ormierten Ebene wird gemaft einer be- 
vorzugten Variante der Erfindung dadurch ermittelt, indem ei- 
ne Konturf ortsetzung in Richtung der Winkelkoordinate (<|>) mit 
einer Optimierung nach minimalen Anderungen erfolgt. 

5 

Es wir also eine Transformation nach Polar koordinaten beziig- 
lich des Volumenschwerpunktes und der ersten Gerade angewen- 
det. Damit ergibt sich ein "Konturf eld" K A (x\<|>) m it festge- 
legtem <|>, das 360° abdeckt. Entsprechend den Polarkoordinaten 

10 wird das radiale x' im Folgenden als r bezeichnet (r-x*). Mit 
relativ geringem Aufwand bleibt nur das "eindimensionale" 
^^^\Problem, den Radius r zu jedem Winkel § als einzige Grofie zu 
bestimmen. Dass dies tatsachlich ein eher „eindimensionales x " 
Problem ist, kann man sich dadurch veranschaulichen, dass man 

15 diese Funktion kartesisch auftragt, z.B. mit (f) in der Abszis- 
se und r in der Ordinate. Fur r ergeben sich nun einige Rand- 
bedingungen, die man sich zunutze machen kann, urn eine quali- 
tativ gute Kontur zu finden: 

20 - Die Kontur muss fur fa und fa die gleiche r - Koordinate 
haben (Geschlossenheit der Kontur) . 

- Die Kontur ist kontinuierlich und hat uberwiegend vertika- 
len Verlauf. Der r-Versatz von Zeile zu Zeile ist relativ 
gering. 

'4fct?5 - Ein typisches "Kantenprof il " im Helligkeit sprof il lasst 
sich durch Zeilenvergleich finden. 

- Moglichkeiten der Kompensation von Signalrauschen im Zei- 
lenvergleich . 



30 Fur die Transformation nach Polarkoordinaten sollte das Win- 
kelinkrement bevorzugt so gewahlt werden, dass benachbarte 
Punkte ((|)i, r) und (c|> i+1 , r) im Bereich der Kontur etwa den 
Abstand einer Voxelkantenlange haben. Anstelle der Transfor- 
mation nach Polarkoordinaten kann auch eine Transformation 

35 beziiglich einer "Modellkontur " erfolgen, wenn eine Ahnlich- 
keit zu einem solchen Modell zu erwarten ist. Die Abszisse 
ist dann durch die f ortlauf enden Indizes ii bis i n der Kon- 



200211387 



9 

turpunkte des Modells gekennzeichnet , auf der Ordinate sind 
die senkrechten Abstande zur Modellkontur aufgetragen. Es 
handelt sich hier wie bei der Transformation nach Polkoordi- 
naten um ein Resampling. 

5 

In der Ergebnismatrix des Resamplings konnen nun verschiedene 
Merkmale fur eine Objektkante bestimmt bzw. durch entspre- 
chende Signalverarbeitung hervorgehoben, sowie in einen Zu- 
sammenhang gebracht werden. Dabei kann z.B. ein Suchverf ahren 
10 angewandt werden, das alle sinnvolle Pf adver zweigungen analy- 

siert und zusammenf ugt . Beispiel einer Realisierung ist die 
!^"\dynamische Optimierung nach R. Bellman, "Dynamic programming 
and stochastic control processes", Information and Control, 
1(3), Seiten 228-239, September 1958. Das Ergebnis ist eine 
15 fortlaufende Kontur, z.B. in der Polarkoordinatendarstellung 
von fyi bis (|) n (bzw. fur i x bis i n ) . Diese Kontur kann in die 
ursprunglichen kartesischen Koordinaten zuriicktransf ormiert 
werden. Die dynamische Optimierung kann mehrere Mer kmalswerte 
(Filterantworten) zu einer geeigneten Kostenf unktion kombi- 
20 nieren und enthalt nach der Vorwarts-Akkumulation von Kosten 
als wesentliches Element ein Backtracking. Damit ist eine E- 
bene des ersten Sets von Ebenen konturiert. 

Entsprechend werden die restlichen Ebenen des ersten Sets von 
\|P'5 Ebenen analysiert. Dabei konnen bereits bestimmte Konturen 
vorher analysierter Ebenen verwendet werden, um auch in der 
durch a beschriebenen Rotationsrichtung die Konsistenz der 
dreidimensionalen Oberf lachenkontur des dreidimensionalen Ab- 
bildes zu gewahrleisten ( verwendbares Vorabwissen) . Ebenso 
30 konnen alle Ebenen des ersten Sets von Ebenen unabhangig ana- 
lysiert werden. Mit einer Oder mehreren nachgeschalteten Be- 
arbeitungsstuf en, mit denen die einzelnen Konturen der Ebenen 
des ersten Sets von Ebene verglichen werden, kann das erhal- 
tene Oberf lachengitter verfeinert werden. 



35 



Am Ende der Bearbeitung aller Ebenen des ersten Sets von Ebe- 
nen liegt die Oberf lachenkontur des dreidimensionalen Abbil- 
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des zunachst in Koordinaten (r, a, <j>) vor. In diesem "Quader" 
ist eine Bearbeitung des 3D-Zusammenhanges moglich, z.B. als 
Filterung der r-Werte mit einem 2D-0perator. 

Die Riicktransf ormation in die ursprtinglichen kartesischen Ko- 
ordinaten ergibt das Oberf lachengitter des dreidimensionalen 
Abbildes. Das Oberf lachengitter ist vergleichbar mit einem 
netzahnlichen Kafig. Das Oberf lachengitter kann nun in den 
Original-Koordinaten verbessert und verfeinert werden. Ab- 
schliefiend konnen die Voxel innerhalb des Oberf lachengitters 
als zum dreidimensionalen Abbild zugehorig markiert werden. 

\ 

i 

Das an einer spharischen Geometrie orientierte erf indungsge- 
mafie Verfahren hat z.B. den Vorteil, dass es keine Vorzugs- 
richtung besitzt und dass, statistisch betrachtet, praktisch 
auftretende Oberf lachenverlaufe durchschnittlich besser be- 
handeln werden als bei einem Verfahren, das sich am kartesi- 
schen Gitter orientiert, da eine solche Oberflache besser pa- 
rallel zu einer geschlossenen 3D-Kontur liegt. Weil das er- 
f indungsgemafte Verfahren von sich aus geometrisch relativ gut 
an die Problemstellung angepasst ist, kann die Anzahl der E- 
benen des ersten Sets von Ebenen um ein Mehrfaches kleiner 
sein, als wenn im Bereich d.er Oberflache im Voxelabstand ab- 
getastet werden wiirde; d.h. die "Net zlinien" konnen bei- 
spielsweise zehn oder mehr Voxelkanten Abstand voneinander 
haben, um die Gesamtkontur immer noch hinreichend genau zu 
beschreiben . 

Insbesondere fur eine verbesserte Ermittlung der Oberflachen- 
konturierung des dreidimensionalen Abbildes oder wenn unguns- 
tige Ausgangsdaten vorliegen, werden gemaft einer Ausfuhrungs- 
form der Erfindung zusatzlich folgende Verf ahrensschritte 
durchgef uhrt : 

- Legen eines zweiten Sets von Ebenen durch das dreidimensio- 
nale Abbild, wobei sich die Ebenen des zweiten Sets von E- 
benen in einer zweiten Gerade schneiden, 
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- Ermitteln der Konturen, die in jedem der Ebenen des zweiten 
Satzes von Ebenen abgebildet und der Oberf lachenkontur des 
dreidimensionalen Abbildes zugeordnet sind, und 

- zusammen mit den in jeder Ebene des ersten Sets von Ebenen 
ermittelten Konturen, Zusammenf iigen der in jeder Ebene des 
zweiten Sets von Ebenen ermittelten Konturen zu dem die O- 
berf lachenkontur des dreidimensionalen Abbildes zugeordne- 
ten Oberf lachengitters . 



10 



Durch Re-Extraktion orientiert an den mittels der Ebenen des 
A ersten Sets von Ebenen ermittelten und dem Oberf lachengitter 
zugeordneten Netzlinien (Scanlinien rechtwinklig zueinander 
verlaufend mit Orientierung von innen nach aufien) gelingt 
15 auch die Konturierung weniger konvexer, unregelmaftiger drei- 
dimensionaler Abbilder. 

Ein Gewinn an Zuverlassigkeit und Genauigkeit kann erzielt 
werden, wenn, wie es nach einer Ausf iihrungsf orm der Erfindung 
20 vorgesehen ist, die erste Gerade rechtwinklig zur zweiten Ge- 
rade ausgerichtet ist und insbesondere die erste Gerade und 
die zweite Gerade sich im Volumenschwerpunkt des dreidimensi- 
onalen Abbildes schneiden. 

'£03 Ein Nutzen der genauen Oberf lachenkontur ierung besteht z.B. 

in der Moglichkeit, das dreidimensionale Abbild genau fur ei- 
ne funktionelle Beurteilung zu bestimmen. Fur therapeutische 
Oder chirurgische Maftnahmen ist die genaue Form in der Bezie- 
hung zur Nachbarschaf t (z.B. Gefafie) und auch die Formveran- 
30 derung im Laufe der Zeit sehr wichtig. Weitere Anwendungen 

sind in der plastischen Chirurgie oder in der Prothetik mog- 
lich im Sinne der Erfassung der zu ersetzenden Objekte oder 
von Kavitaten fur ein Implantat. 



35 



Eine vereinfachte Berechnung der benotigten Rechenoperationen 
soil in ungefahr aufzeigen, welchen Vorteil das erf indungema- 
fte Verfahren hat. Es wird angenommen, dass der Volumendaten- 
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satz in Form von 50 Computertomographie-Schichten der Matrix 
512x512 vorliegt. Das dreidimensionale Abbild ist in einem 
Teilquader von 40x200x200 Voxeln enthalten. 

5 Bei einer Oberf lachenkonturierung des dreidimensionalen Ab- 
bildes gemaft dem Stand der Technik wird ausgehend von einem 
Startpunkt auf der Oberf lache deren Fortsetzung gesucht. Die 
Oberf lache des gesamten Teilquaders betragt 112000 "Pixel". 
Fur die Gradientenbildung wird ein Faktor 11 bezuglich der 

10 Oberf lache angesetzt (Koppelfeld 3 3 — > "26") und mit einem 
Faktor "1" fur die 27 uberdeckten Voxel angesetzt. Dies ent- 
^JpT spricht effektiv funf Oberf lachen-Suchstreifen innerhalb und 
aufierhalb der Oberf lache. Von jedem Punkt ausgehend werden in 
einer Vor zugsrichtung und seitlich im Prinzip 17 Richtungen 

15 untersucht (anfangs) , die sich dann innerhalb einer Flache 

reduzieren. Sinnvollerweise begrenzt man die Suche empirisch, 
wie z.B. in Buckley K., Vaddiraju A., Perry R. , "A New Pru- 
ning/Merging Algorithm for MHT Mult itarget Tracking", Radar- 
2000, Mai 2000 beschrieben. Ferner wird angenommen, dass es 

20 geniigt, durchschnittlich sechs Alternativen mit einer durch- 
schnittlichen Tiefe von drei zu untersuchen, urn eine sinnvol- 
le Fortsetzung von einem Irrweg zu unterscheiden bzw. urn auf 
einem schon klassif izierten Punkt zu landen (Mult iplikator 
6 3 =216) . 

# 

Bei einer Oberf lachenkonturierung gemafi dem eben beschriebe- 
nen Verfahren und unter nicht Berucksichtigung des sogenannte 
"Backtrackings", benotigt man in etwa 11x112000x26=33264000 
Gradientenbildungen und fur die Suche nach maximaler Gradien- 
30 tensumme in etwa 216x112000=24192000 Rechenoperationen . Somit 
ergeben sich in etwa 57000000 Rechenoperationen. 

Fur die ungefahre Berechnung benotigter Rechenoperationen fur 
die Oberf lachenkonturierung desselben dreidimensionalen Ab- 
35 bildes gemaft dem erf indungsgemaften Verfahrens wird angenom- 
men, dass das erste Set von Ebenen 20 Ebenen umfasst. Nach 
der Transformation in Polar koordinaten werden die Bildinfor- 
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mationen der transf ormierten Ebenen jeweils im 2°-Raster ab- 
getastet, wobei die entsprechenden CT-Bilder im ublichen Ras- 
ter 512x512 rekonstruiert werden: 180x256=46000. Die Gradien- 
tenbildung und Kontursuche werden zunachst in den Ebenen 
5 durchgef uhrt . Aufierdem wird eine rauschunterdruckende Fal- 
tungsoperation mit einem Koppelfeld von sieben durchgef uhrt . 
Dabei wird ein typisches Dichteprofil im Bereich des Objekt- 
randes mit sieben Stutzstellen erzeugt und uber den gesamten 
Radius verschoben. in jeder Verschiebeposition wird die Kova- 

10 rianz als Iterationswert errechnet. Bei der Suche der Kontu- 
ren ist jeweils ein "Sprung" urn 1 Pixel von Zeile zu Zeile 
vorgesehen, es gibt also drei Alternativen . Zur Nachbearbei- 
tung des Zusammenhangs zwischen den Ebenen wird ein 3x3- 
Koppelfeld verwendet . Es ergeben sich somit in etwa 

15 20x46000=920000 Operationen fur die Gradientenbildung, und 
fur die Suche nach maximaler Gradientensumme ca. 
3x4 6000=138000 Rechenoperationen. Fur das Backtracking, das 
hier einer Uberprufung des Ebenen-Zusammenhangs entspricht, 
ergeben sich ungefahr 20x180x9=32400 Rechenoperation . Somit 

20 ergeben sich in etwa 1090400 Rechenoperationen. 

Aufgrund des erfindungsgemaften Verfahrens mit den genannten 
Annahmen ergibt sich eine Verringerung des Aufwandes an Re- 
chenoperationen urn etwa den Faktor 50. 

Ein Ausf uhrungsbeispiel ist exemplarisch in den beigefugten 
schematischen Zeichnungen dargestellt. Es zeigen: 

Figur 1 einen Computertomographen, 

30 

Figur 2 ein Set von Ebenen, deren Ebenen ein mit dem in der 

Figur 1 dargestellten Computertomographen aufgenomme- 
nes dreidimensionales Abbild schneiden, 

35 Figur 3 in kartesischen Koordinaten dargestellte Bildinforma- 
tionen einer der in der Figur 2 dargestellten Ebenen, 
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Figur 4 in Polarkoordinaten dargestellte Bildinf ormationen 
der in der Figur 3 dargestellten Ebene und 

Figur 5 eine Veranschaulichung eines Ermittelns einer in der 
Figur 4 dargestellten Kontur. 

Die Figur 1 zeigt schematisch einen Computertomographen mit 
einer Rontgenstrahlenquelle 1, von dem ein pyramidenf ormiges 
Rontgenstrahlenbundel 2, dessen Randstrahlen in der Figur 1 
strichpunktiert dargestellt sind, ausgeht, das ein Untersu- 
chungsobjekt , beispielsweise einen Patienten 3, durchsetzt 
und auf einen Strahlungsdetektor 4 trifft. Die Rontgenstrah- 
lenquelle 1 und der Strahlungsdetektor 4 sind im Falle des 
vorliegenden Ausf iAhrungsbeispieles an einer ringf ormigen 
Gantry 5 einander gegenuberliegend angeordnet. Die Gantry 5 
ist bezuglich einer Systemachse 6, welche durch den Mittel- 
punkt der ringformigen Gantry 5 verlauft, an einer in der Fi- 
gur 1 nicht dargestellten Halterungsvorrichtung drehbar ge- 
lagert (vgl. Pfeil a). 

Der Patient 3 liegt im Falle des vorliegenden Ausf iihrungsbei- 
spieles auf einem fur Rontgenstrahlung transparenten Tisch 7 , 
welcher mittels einer in der Figur 1 ebenfalls nicht darge- 
stellten Tragevorrichtung langs der Systemachse 6 verschieb- 
bar gelagert ist (vgl. Pfeil b) . 

Die Rontgenstrahlenquelle 1 und der Strahlungsdetektor 4 bil- 
den somit ein Messsystem, das bezuglich der Systemachse 6 
drehbar und entlang der Systemachse 6 relativ zum Patienten 3 
verschiebbar ist, so dass der Patient 3 unter verschiedenerK 
Projektionswinkeln und verschiedenen Positionen bezuglich der 
Systemachse 6 durchstrahlt werden kann. Aus den dabei auftre- 
tenden Ausgangssignalen des Strahlungsdetektors 4 bildet ein 
Datenerf assungssystem 9 Messwerte, die einem Rechner 11 zuge- 
fuhrt werden, der mittels dem Fachmann bekannten Verfahren 
ein Bild des Patienten 3 berechnet, das wiederum auf einem 
mit dem Rechner 11 verbundenen Monitor 12 wiedergegeben wer- 
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den kann . Im Falle des vorliegenden Ausf iihrungsbeispiels 1st 
das Datenerf assungssystem 9 mit einer elektrischen Leitung 8, 
die in nicht dargestellter Weise beispielsweise ein Schleif- 
ringsystem oder eine drahtlose Ubertragungsstrecke enthalt, 
5 mit dem Strahlungsdetektor 4 und mit einer elektrischen Lei- 
tung 14 mit dem Rechner 11 verbunden. 

Das in der Figur 1 gezeigte Computertomographiegerat kann so- 
wohl zur Sequenzabtastung als auch zur Spiralabtastung einge- 
10 setzt werden. 



Bei der Sequenzabtastung erfolgt eine schichtweise Abtastung 
des Patienten 3. Dabei wird die Rontgenstrahlenquelle 1 und 
der Strahlungsdetektor 4 bezuglich der Systemachse 6 urn den 
15 Patienten 3 gedreht und das die Rontgenstrahlenquelle 1 und 
den Strahlungsdetektor 4 umfassende Messsystem nimmt eine 
Vielzahl von Projektionen auf, um eine zweidimensionale 
Schicht des Patienten 3 abzutasten. Aus den dabei gewonnen 
Messwerten wird ein die abgetastete Schicht darstellendes 
20 Schnittbild rekonstruiert . Zwischen der Abtastung aufeinan- 

derfolgender Schichten wird der Patient 3 jeweils entlang der 
Systemachse 6 bewegt. Dieser Vorgang wiederholt sich so lan- 
ge, bis alle interessierenden Schichten erfasst sind. 

15 Wahrend der Spiralabtastung dreht sich das die Rontgenstrah- 
lenquelle 1 und den Strahlungsdetektor 4 umfassende Messsys- 
tem bezuglich der Systemachse 6 und der Tisch 7 bewegt sich 
kontinuierlich in Richtung des Pfeils b, d.h. das die Ront- 
genstrahlenquelle 1 und den Strahlungsdetektor 4 umfassende 
30 Messsystem bewegt sich relativ zum Patienten 3 kontinuierlich 
auf einer Spiralbahn c, so lange, bis der interessierende Be- 
reich des Patienten 3 vollstandig erfasst ist. Dabei wird ein 
Volumendatensat z generiert, der im Falle des vorliegenden 
Ausf iihrungsbeispiels nach dem in der Medizintechnik ublichen 
35 DICOM-Standard kodiert ist. 
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Im Falle des vorliegenden Ausf uhrungsbeispiels soil die Ober- 
flache eines mit dem Volumendatensat z aufgenommenen dreidi- 
mensionalen Abbildes 21 eines Objektes, das Teil des interes- 
sierende Bereichs des Patienten 3 ist, konturiert werden. Im 
Falle des vorliegenden Ausf uhrungsbeispiels handelt es sich 
bei dem Objekt urn eine Gewebsveranderung des Patienten 3. Das 
dreidimensionale Abbild 21 ist in der Figur 2 schematisch 
dargestellt und kann mit dem Monitor 12 betrachtet werden. 

Urn die Oberflache des dreidimensionalen Abbildes 21 zu kontu- 
rieren, wird ein erstes Set Ebenen, das im Falle des vorlie- 
genden Ausf uhrungsbeispiels 20 Ebenen Ei bis E 2 o umfasst, 
durch das dreidimensionale Abbild 21 gelegt. Der Obersicht 
halber sind in der Figur 2 nur die Ebene Ei bis E 5 darge- 
stellt. Die Ebenen E\ bis E 2 o schneiden sich in einer ersten 
Geraden G, die im Falle des vorliegenden Ausf uhrungsbeispiels 
durch den Volumenschwerpunkt S des dreidimensionalen Abbildes 
21 verlauft. Der Volumenschwerpunkt wird entweder manuell von 
einer in der Figur 1 nicht dargestellten Bedienperson mittels 
einer in der Figur 2 nicht dargestellten, in das Sichtfeld 
des Monitors 12 einblendbare Marke, die mit einer mit dem 
Rechner 11 verbundenen Rechnermaus 13 bewegt werden kann, im 
dreidimensionalen Abbild 21 ungefahr markiert oder der Rech- 
ner 11 ist im Falle des vorliegenden Ausf uhrungsbeispiels 
derart konf iguriert , dass er automatisch den ungefahren Volu- 
menschwerpunkt S des dreidimensionalen Abbildes 21 ermittelt 
und anschliefiend automatisch die Ebenen Ei bis E 2 o bestimmt. 
Aufterdem kann die Bedienperson im Falle des vorliegenden Aus- 
f uhrungsbeispiels den Volumendatensat z mit dem Rechner 11 
beispielsweise einer Medianf ilterung unterziehen. Der Volu- 
mendatensatz wird dadurch unter Kantenerhaltung stark geglat- 
tet, urn die automatische Ermittlung des Volumenschwerpunkts 
zu verbessern. 

Im Falle des vorliegenden Ausf uhrungsbeispiels sind die Ebe- 
nen Ei bis E 2 o urn jeweils einen Winkel Aa=9° gegeneinander 
verdreht und die erste Gerade G verlauft parallel zur x-Achse 
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des fur den Volumendatensat z bereits definierten DICOM- 
Koordinatensystems . Umfasst das dreidimensionalen Abbild 21 
beispielsweise ein Viertel der Gesamtgrofie des im Falle des 
vorliegenden Ausf uhrungsbeispiels isotropen Volumendatensat- 
5 zes (512er Kubus) , so entspricht der maximale Punktabstand an 
der Oberflache des dreidimensionalen Abbildes 21 in etwa 20 
Pixel . 



In jedem der Ebenen Ei bis E 2 o sind Schnitte des dreidimensio- 
10 nalen Abbildes 21 abgebildet. Ein solcher Schnitt 31, der im 

Falle des vorliegenden Ausf uhrungsbeispiels durch die Ebene 
K Ei festgelegt ist, ist schematisch in der Figur 3 gezeigt. 

Aufterdem sind die Bildinhalte der einzelnen Ebenen Ei bis E 2 o 

in kartesischen Koordinaten beschrieben. Die Bildinhalte der 
15 Ebene Ei sind im Falle des vorliegenden Ausf uhrungsbeispiels 

mit den kartesischen Koordinaten x x> und y^' beschrieben. 

Urn die Oberf lachenkontur des dreidimensionalen Abbildes 21 zu 
erhalten, werden anschliefiend die in den Ebenen Ei bis E 2 o 
20 dargestellten Schnitte, also auch der Schnitt 31 der Ebene Ei 
analysiert, indem die Konturen der einzelnen Schnitte ermit- 
telt werden. Die Kontur des Schnittes 31 ist in der Figur 3 
mit dem Bezugszeichen 32 versehen. 

Fur eine effiziente Ermittlung der Konturen werden im Falle 
des vorliegenden Ausf uhrungsbeispiels die Bildinf ormationen 
jeder Ebene Ei bis E 2 o radial neu angeordnet, indem sie bezug- 
lich des Volumenschwerpunktes S nach Polarkoordinaten (r, 
transf ormiert werden. Das Ergebnis einer solchen Koordinaten- 
30 transformation ist in der Figur 4 beispielhaft dargestellt. 
Es handelt sich dabei urn die Transformation der Bildinhalte 
und insbesondere des in der Figur 3 dargestellten Schnittes 
31 der Ebene Ei. Die transf ormierte Kontur des Schnittes 31 
ist in der Figur 4 mit dem Bezugszeichen 32 1 versehen. 

35 

Anschliefiend wird die transf ormierte Kontur 32 1 ermittelt. Es 
wird also ein Begrenzungspf ad gesucht, der der in Polarkoor- 
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dinaten transf ormierten Kontur 32 1 entspricht. Dazu erfolgt 
zunachst eine richtungsbetonte Filterung, die sogar eindimen- 
sional in r-Richtung sein kann und die sowohl kantensensitiv 
als auch rauschunterdruckend wirkt. Uber die Filterungen 
5 lauft nun nach dem Prinzip der dynamischen Optimierung, wie 
z.B. in R. Bellman, "Dynamic programming and stochastic 
control processes", Information and Control, 1(3), Seiten 
228-239, September 1958 beschrieben, eine eindimensionale Su- 
che nach dem Begrenzungspf ad, also nach der in Polarkoordina- 

10 ten transf ormierten Kontur 32 ? . Es wird also fur die Winkel- 
koordinaten § zwischen 0° und 360° entsprechenden r-Koordinaten 
^ gesucht, wobei als Randbedingung gilt, dass die r-Koordinate 
fur <|)=0 o gleich der r-Koordinate flir <|>=360 o sein soil. Diese 
r-Koordinaten werden im Falle des vorliegenden Ausfiihrungs- 

15 beispiels ermittelt, indem nach einer relativ starken Glat- 
tung die ersten und letzten funf Zeilen iiberlagert und an- 
schlieftend die Kantenposition in den geglatteten Zeilendaten 
ermittelt werden. Diese Suche ist exemplarisch in der Figur 5 
gezeigt, wobei die transf ormierte Kontur wieder mit dem Be- 

20 zugszeichen 32 1 versehen ist. In der Figur 5 gezeigte Linien 
40 sind dunner als die transf ormierte Kontur 32' gezeichnet 
und reprasentieren wahrend der Ermittlung gefundene Pfadal- 
ternativen, die jedoch wieder verworfen wurden. 

I 

|25 Ahnlich werden die in Polarkoordinaten transf ormierte Kontu- 
ren der restlichen Ebenen E 2 bis E 2 o ermittelt und anschlie- 
ftend in Zusammenhang betrachtet. Darunter wird verstanden, 
dass untersucht wird, ob die Kontinuitat der Helligkeit bzw. 
des in dieser Richtung relevanten Gradienten plausibel sind. 

30 

Anschlieliend werden die gefundenen in Polarkoordinaten trans- 
f ormierten Konturen jeder Ebene E x bis E 2 o zuriick in kartesi- 
sche Koordinaten transf ormiert und danach zu einem Oberfla- 
chengitter zusammengef iigt . Das Oberf lachengitter entspricht 
35 dabei der Oberf lachenkontur des dreidimensionalen Abbildes 
21. 
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Im Falle des vorliegenden Ausf uhrungsbeispiels wird nun das 
Oberf lachengitter in den ursprunglichen dreidimensionalen 
kartesischen Koordinaten verbessert und verfeinert. Anhand 
des Oberf lachengitters wird die zu ermittelnde Oberflachen- 
5 kontur des dreidimensionalen Abbildes 21 bis auf einfachen 
Pixelabstand interpoliert . Abschliefiend werden die Voxel in- 
nerhalb des Oberf lachengitters als zum dreidimensionalen Ab- 
bild 21 zugehorig markiert. 

10 Durch das Legen der Ebenen Ei bis E 2 o durch das dreidimensio- 
nale Abbild 21 und durch die anschlieftende Koordinatent rans- 
formation nach Polarkoordinaten wird der Volumendatensat z als 
dreidimensionaler Datenblock B'(r, § r a) dargestellt. Das Er- 
mitteln der transf ormierten Konturen, wie z.B. der transfor- 

15 mierten Kontur 31 f , ist fur nahezu konvexe Objekte weniger 
rechenauf wandig . Es wird eine zusammenhangende "Gebirgsfla- 
che" uber <j> und a gesucht. Die "Gebirgsf lache" weist lokal 
relativ geringe Schwankungen in r-Richtung auf. 

20 Auch die Anzahl der fur die Ermittlung der Konturen analy- 

sierten Voxel wird gegenuber einer dreidimensionalen Oberfla- 
chenkonturierung gemafi dem Stand der Technik verringert. Um- 
fasst der Volumendatensat z z.B. 50 CT-Schichten mit einer un- 
m gefahren Ausdehnung von 128 Pixeln des dreidimensionalen Ab- 

(j(p25 bildes f dessen Oberf lachenkontur ermittelt werden soil, also 
einem etwa 50x256x256 groften Bearbeitungsf eld, werden in etwa 
eine Million Voxel (20 Ebenen, Polar koordinatenraster 2°, r max 
= 256) analysiert. Gemafi dem Stand der Technik mussten dage- 
gen in etwa drei Millionen Voxel analysiert werden. 

30 

Urn das verbessertes Oberf lachengitter zu erhalten, ist es im 
Falle des vorliegenden Ausf uhrungsbeispiels vorgesehen, dass 
zusatzlich zu den Ebene Ei bis E 2 o des ersten Sets von Ebenen 
wahlweise wenigstens ein weiteres, in den Figuren der Uber- 
35 sicht halber nicht dargestelltes zweites Set von Ebenen ge- 
legt werden kann. Die Ebenen des zweiten Sets von Ebenen 
schneiden sich in einer zweiten Gerade, die ebenfalls durch 



200211387 



20 

den Volumenschwerpunkt S verlauft und rechtwinklig zur ersten 
Gerade G ausgerichtet ist. Die mit diesen Ebenen dargestell- 
ten Schnitte werden wieder in Polarkoordinaten bezuglich des 
Volumenschwerpunktes S transf ormiert . Anschlieftend werden 
wieder die transf ormierten Konturen, die der Oberf lachenkon- 
tur des dreidimensionalen Abbildes 21 zugeordnet sind, analog 
wie obenstehend ermittelt. Die ermittelten transf ormierten 
Konturen werden anschlieftend wieder zuruck transf ormiert und 
danach zu dem Oberf lachengitter zusammengef ugt . 

Ein Vorteil ist auch, dass durch die Koordinatentransf ormati- 
on in Polarkoordinaten das Ermitteln der transf ormierten Kon- 
turen in Drehrichtung a (azimutal) und in Richtungen inner- 
halb der Ebenen (<j), orbital) im Wesentlichen als zweidimensi- 
onale Filterung, namlich als Glattung von r(ty,a), anstelle in 
Originalkoordinaten dreidimensional durchgefuhrt werden kann. 

Der Volumendatensatz, mit dem das dreidimensionale Abbild 21 
hergestellt wird, wird im Falle des vorliegenden Ausfuhrungs- 
beispiels mit einem Computertomographen hergestellt. Der Vo- 
lumendatensatz kann jedoch auch mit anderen medizintechni- 
schen Geraten, wie insbesondere mit einem Rontgengerat , einem 
Magnetresonanzgerat, einem PET, einem Ultraschallgerat oder 
mittels Nuklearmedizin hergestellt werden. Auch bildgebende 
Gerate oder technische Einrichtungen aufierhalb der Medizin- 
technik sind verwendbar. 

Das dreidimensionale Abbild 21 ist ein Abbild eines Teiles 
des Patienten 3. Mittels des erf indungsgemaften Verfahrens 
kann aber auch die Oberflache eines nichtmenschlichen dreidi- 
mensionalen Abbildes konturiert werden. 

Die erste Gerade G muss auch nicht notwendigerweise durch den 
Volumenschwerpunkt S des dreidimensionalen Abbildes 21 ver- 
laufen. Die Ebenen Ei bis E 2 o mussen sich auch nicht notwendi- 
gerweise in der Geraden G schneiden. 
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Patent anspruche 

1. Verfahren zur Oberf lachenkonturierung eines dreidimensio- 
nalen Abbildes (21) eines Objekts, aufweisend folgende Ver- 
f ahrensschritte : 

- Legen eines ersten Sets von Ebenen (Ei-E 2 o) durch das drei- 
dimensionale Abbild (21) , 

- Ermitteln der Konturen (32) , die in jeder der Ebenen (Ei~ 
E 2 o) des ersten Satzes von Ebenen abgebildet und der 
Oberf lachenkontur des dreidimensionalen Abbildes (21) 
zugeordnet sind, und 

- Zusammenf iigen der in jeder Ebene (Ei-E 2 o) des ersten Sets 
von Ebenen ermittelten Konturen (32) zu einem die Oberfla- 
chenkontur des dreidimensionalen Abbildes (21) zugeordneten 
Oberf lachengitters . 

2. Verfahren nach Anspruch 1, bei dem sich die Ebenen (Ei-E 2 o) 
des ersten Sets von Ebenen in einer ersten Gerade (G) schnei- 
den . 

3. Verfahren nach Anspruch 2, bei dem die erste Gerade (G) 
wenigstens im Wesentlichen durch den Volumenschwerpunkt (S) 
des dreidimensionalen Abbildes (21) verlauft. 

4. Verfahren nach Anspruch 3, bei dem zum Ermitteln des Volu- 
menschwerpunkts (S) des dreidimensionalen Abbildes (21) ein 
Volumendatensatz, mittels dessen das dreidimensionale Abbild 
(21) gespeichert ist, einer Filterung unterzogen wird. 

5. Verfahren nach Anspruch 3 oder 4, bei dem der Volumen- 
schwerpunkt (S) des dreidimensionalen Abbildes (21) automa- 
tisch bestimmt wird. 
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6. Verfahren nach einem der Anspruche 1 bis 5, bei dem das 
dreidimensionale Objekt (21) ein Teil eines Lebewesens (3) 
ist und das dreidimensionale Abbild (21) mittels eines medi- 
zintechnischen Gerates aufgenommen wird. 

5 

7. Verfahren nach Anspruch 6, bei dem der Volumenschwerpunkt 
(S) des dreidimensionalen Abbildes (21) wahrend einer naviga- 
tionsgef uhrten Behandlung des Lebewesens (3) ermittelt wird. 

10 8. Verfahren nach einem der Anspruche 3 bis 7, bei dem die 
Bildinhalte jeder Ebene (Ei-E 2 o) des ersten Sets von Ebenen 

^ mit kartesischen Koordinaten (x", y" ) beschrieben sind und 
fur das Ermitteln der Konturen (32) jeder einzelnen Ebene 
(Ei-E 2 o) folgende Verf ahrensschritte fur jede Ebene (Ei~E 2 o) 

15 ausgefiihrt werden: 

- Durchfuhren einer Koordinatentransf ormat ion nach Polarkoor- 
dinaten (r,<|>), beziiglich des Volumenschwerpunktes (S) des 
dreidimensionalen Abbildes und damit Abrollen der Kontur 

20 (32) und 

- Ermitteln der Kontur (32^) in der transf ormierten Ebene. 

9. Verfahren nach Anspruch 8, bei dem die Kontur (32 x ) in der 
25 transf ormierten Ebene dadurch ermittelt wird, indem eine Kon- 

turf ortset zung in Richtung der Winkelkoordinate (())) mit einer 
Optimierung nach minimalen Anderungen erfolgt. 

10. Verfahren nach Anspruch 8 oder 9, bei dem die Kontur 

30 (32M in der transf ormierten Ebene mittels dynamischer Opti- 
mierung verbessert wird. 

11. Verfahren nach einem der Anspruche 2 bis 10, bei dem zu- 
satzlich folgende Verf ahrensschritte durchgefiihrt werden: 

35 
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- Legen eines zweiten Sets von Ebenen durch das dreidimensio- 
nale Abbild (21) , wobei sich die Ebenen des zweiten Sets 
von Ebenen in einer zweiten Gerade schneiden, 

5 - Ermitteln der Konturen, die in jedem der Ebenen des zweiten 
Satzes von Ebenen abgebildet und der Oberf lachenkontur des 
dreidimensionalen Abbildes zugeordnet sind, und 

- zusammen mit den in jeder Ebene (Ei-E 2 o) des ersten Sets von 
10 Ebenen ermittelten Konturen, Zusammenf ugen der in jeder E- 



15 12. Verfahren nach Anspruch 11, bei dem die erste Gerade (G) 
rechtwinklig zur zweiten Gerade ausgerichtet ist. 

13. Verfahren nach Anspruch 12, bei dem sich die erste Gerade 
(G) und die zweite Gerade im Volumenschwerpunkt (S) des drei- 
20 dimensionalen Abbildes (21) schneiden. 



bene des zweiten Sets von Ebenen ermittelten Konturen zu 
dem die Oberf lachenkontur des dreidimensionalen Abbildes 
(21) zugeordneten Oberf lachengitters . 
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Zusammenf as sung 

Verfahren zur Oberflachenkonturierung eines dreidimensionalen 
Abbildes 

Die Erfindung betrifft ein Verfahren zur Oberflachenkonturie- 
rung eines dreidimensionalen Abbildes (21) eines Objekts. Zu- 
erst wird ein erstes Set von Ebenen (Ei-E 2 o) durch das dreidi- 
mensionale Abbild (21) gelegt und die Konturen (32), die in 
jeder der Ebenen (Ei-E 20 ) des ersten Satzes von Ebenen abge- 
bildet und der Oberf lachenkontur des dreidimensionalen Abbil- 
des (21) zugeordnet sind, ermittelt. Die Konturen (32) werden 
anschliefiend zu einem die Oberf lachenkontur des dreidimensio- 
nalen Abbildes (21) zugeordneten Oberf lachengitters zusammen- 
gelegt . 



FIG 2 
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